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.’n  Unlijuited 

A LIAPUNOV  FUNCTIONAL  FOR  A MATRIX 
NEUTRAL  DIFFERENCE-DIFFERENTIAL  EQUATION  WITH  ONE  DELAY 

V 

Abstract;  For  the  matrix  neutral  difference-differential  equation 
x^(t)  + Ax^(t-l)  = Bx(t)  + Cx(t-T)  we  eonstruet,  a quadratic 

'A*  .4  (■  »..  I', f ■' t- .( 

Liapunov  functional^ which  gives  necessary  and  sufficient  conditions 
for  the  asymptotic  stability  of  the  solutions  of  that  equation. 

We  con-sdrder  a difference  equation  approximation^  of  the  difference- 
differential  equation,  and  for  this  difference  equation  we  construct 

,i  Co.  .J  If 

a Liapunov  function  from  which  we,  obtain«»*the  desired  Liapunov  func- 
tional  by  an  appropriate  limiting  process.  The  Liapunov  functional 
thus  obtained  gives  the  best  possible  estimate  for  the  rates  of 
growth  or  decay  of  the  solutions  of  the  matrix  neutral  difference- 
differential  equation.  The  results  obtained  are  natural  generaliza- 
tions of  previous  results  obtained  for  a matrix  retarded  difference- 
differential  equation  with  one  delay. 

1.  Introduction, 

Consider  the  linear  autonomous  matrix  neutral  difference- 
differential  equation  with  one  delay 


□ □ 


asymptotic  behavior  of  the  solutions  of  equation  (1.1). 

For  a linear  autonomous  retarded  difference-differential 
equation  with  one  delay  this  same  problem  was  considered 
recently  in  [9]  for  the  scalar  case,  and  in  [8]  for  the  matrix 
case.  We  will  show  that  the  procedures  developed  in  [9]  and 
[8]  can  be  extended  for  the  neutral  difference-differential 
equation  (1.1) . 

Several  authors,  notably  Hale  [3,6],  Melvin  [14]  and 
Sendaula  [15],  have  considered  Liapunov  functionals  for  neutral 
functional  differential  equations.  Those  Liapunov  functionals 
only  give  sufficient  conditions  for  the  asymptotic  stability  of 
the  solutions.  The  Liapunov  functional  constructed  in  this  paper 
gives  necessary  and  sufficient  conditions  for  the  asymptotic 
stability  of  the  solutions  of  (1.1);  moreover,  it  provides  the 
best  possible  estimate  for  the  rates  of  growth  or  decay  of  the 
solutions. 

As  in  [9] , we  first  consider  a difference-equation  approxima- 
tion of  the  scalar  equation 

x(t)  + ax(t-'t)  = bx(t)  + cx(t-T),  t > 0.  (1.2) 

For  such  a difference  equation  we  construct,  by  means  of  well- 
)cnown  methods,  a Liapunov  function  that  gives  necessary  and 
sufficient  conditions  for  the  asymptotic  stability  of  the  solutions. 
Taking  appropriate  limits  on  this  Liapunov  function  we  obtain  the 
desired  Liapunov  functional  for  equation  (1.2). 
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The  Liapunov  functional  for  the  matrix  neutral  difference- 
differential  equation  (1.1)  is  then  obtained  as  a generalization 
of  the  one  for  the  scalar  equation.  As  in  [8] , the  structure  of 
the  Liapunov  functional  is  completely  analyzed. 

The  Liapunov  functional  constructed  depends  critically 
on  a matrix  function  which  must  satisfy  a special  functional 
differential  equation.  The  existence,  uniqueness  and  structure 
of  the  solutions  of  this  equation  are  described.  This  special 
equation  is  a generalization  of  the  one  studied  in  [1] . 

2 . The  Neutral  Difference-Differential  Equation, 

Let  be  the  space  of  all  Lebesgue  square 

integrable  functions  defined  on  [a,b]  with  values  on  and 

X n 

denote  by  W2([a,b],^  ) the  space  of  all  absolutely  continuous 
functions  which  have  the  first  derivative  in  L2  ( [a  ,b]  , . 

With  T > 0 fixed,  consider  the  Hilbert  space  ^ 0]  , 

with  the  inner  product 

= (|>^(0)())2  (0)  + 1°  <t>i(9)<i)2  (0)d0  + (0  ) 4)  2 (6  ) de 

and  the  induced  norm 

Il'i’lli/  = (0)  + <j)'^(0)(|)  (0)d0  + f c()''^(0)4.'  (G)d0. 

se  J _t  J _x 

Here,  the  superscript  T denotes  the  transpose  of  a matrix. 

Let  x,x:  ->•  ^ ; then  for  t > 0 we  define  the  functions 

^t'*t'  by  x^(0)  = x(t+0),  x^{0)  = x(t+0). 


Consider  the  matrix  neutral  difference-differential  equation 


x(t)  + AxCt-T)  = Bx(t)  + C(t-T),  t > 0, 


(2.1) 


where  A,B,C  are  n x n constant  matrices,  x(t)  is  an  n-vector 
and  T > 0.  Let 


= ♦ 


(2.2) 


be  a given  initial  condition,  with  <p  G 

A solution  of  tha  initial  value  problem  (2.1)-(2.2)  is  a 


function 

x: 

[-T,«.) 

(a) 

^t 

e 

(b) 

X 

satisf 

on 

[0,~) 

(c) 

^0 

= <1>. 

,n 


It  is  ]<nown  that  the  initial  value  problem  (2.1) -(2. 2)  has 
a unique  solution  which  depends  continuously  on  the  initial  data 
in  the  norm  of  A proof  of  this  result  is  given  by  Melvin 

[13] , who  used  the  norm 

1/2 


l<l>  I I = U (-t)  I + 


"^(9)4)'  (0)  l^de 


-X 


which  is  easily  seen  to  be  equivalent  to  the  norm  defined  on  our 
Hilbert  space. 

We  consider  the  solution  operator  T(t):  ^ , defined 


by  T(t)<j>  = x^,  <|)  ^ t > 0.  The  family  T(t),  t > 0,  is  a 
CQ-semigroup.  The  infinitesimal  generator  -Q^  of  T(t)  is 
given  by  = <P ' and  it  has  a domain  , dense  in 

defined  by 

= {<}>  e C ^ and  <t>’  (0)  + A<t>  ' (-T)  = B<t>  (0)  + C<P  (-^)  }, 


[3]  . 


Let  o {jj/)  denote  the  spectrum  of  i.e., 


o(q^)  = {X  I det [X (I+Ae  - B - Ce  = 0}. 


(2.3) 


Then,  [3],  there  exists  a constant  Y such  that  Re(X)  < y for  all 
X G a (ja^)  . Also,  for  every  e > 0 there  exists  a constant 
K > 1 such  that 


I |T(t) 


(Y+e) t 


(2.4) 


For  our  purposes,  it  is  convenient  to  consider  a representa- 
tion of  the  solutions  of  (2.1)  which  is  given  for  every  t,u  > 0 
by  the  formula 


Xf+u^O)  = t^(t)  + Y(t-T)A]x^(0)  + 


J Y (u-ct-T)  [Cx^  (a)  - Ax^(a)]doi, 


(2.5) 
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where  the  matrix  Y is  the  solution  of  initial  value  problem 


Y(t-u)  = I + 


t-T 


t-T 


u 


dgY(t-3-'t)A  + 


Y(t-3)Bd6  + 


Y(t-3--^)Cd3, 


u 


(2.6) 


Y(0)  = I,  Y(t-u)  = 0 for  t < u. 


The  integrals  in  (2.6)  are  Lebesgue-Stiel jes  integrals,  Y(t-u) 
as  a function  of  u is  left  continuous  and  Var  Y(t-.)  < «> 
on  [u,t]  for  every  t > u,  [7].  From  (2.6)  we  have 


Y(t-u)  + Y(t-u-T)A  = I + A + 
rt 


[Y(t-3)B  + Y (t-3-T:)C]d3. 


(2.7) 


u 


For  G it  follows  from  (2.5)  that  Y(t-u)  + Y(t-u-'T)A 

is  an  absolutely  continuous  matrix  valued  function.  Then  from  (2.7) 
we  have  that 


^ [Y(t-u)  + Y(t-u-T)A]  = 
au 


-Y(t-u)B  - Y(t-u-T)C,  a.e.. 


or 


^ [Y(v)  + Y(v-x)A]  = Y(v)B  + Y(v-T)C,  a.e.. 


(2.8) 
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3.  A Liapunov  Function  for  a Difference  Equation  Approximation 
of  the  Scalar  Neutral  Difference-Differential  Equation. 

Our  purpose  in  this  section  is  the  construction  of  a 
Liapunov  function  for  a difference  equation  approximation  of  the 
scalar  neutral  difference-differential  equation 


x(t)  + ax(t-'C)  = bx(t)  + cx(t-T:),  t > 0 . 


(3.1) 


Let  N be  fixed,  and  let  the  intervals  [0,°°),  [-'C,0]  be  sub- 

. . T 

divided  in  subintervals  of  equal  length  — . Denote  the  values  of 

T I T 

the  function  x (9)  at  the  mesh  points  x ^ / by  x^ , 

[ ^ N 

k = 0,1,...,  J = We  thus  obtain  the  difference  equation 


0 ^ , N N,  T 0 , N, 

’'k  ‘‘•’‘k+i  - ' ij  <‘”‘k  ==‘k’ 


J = 1,  . . . ,N, 


an  approximation  to  (3.1).  We  assume,  in  this  section,  that 
c 7^  0,  a ^ 0,  and  rewrite  this  difference  equation  in  the  form 


^k+l  = ^^k' 


(3.2) 


where  denotes  the  (2N+1) -dimensional  vector 


0 1 

0 1 N \ \ 

^k=  CX^/CX,,/..-/CX^r 

N 


.,  - a 


N-1  N 


(3.3) 


and  the  (2N+1)  x (2N+1)  matrix  A is  given  by 


It  is  well-known  how  to  construct  a Liapunov  function  for  the 
difference  equation  (3.2)  so  as  to  obtain  necessary  and  sufficient 
conditions  for  the  asymptotic  stability  of  its  solutions,  [11,12], 
For  this  purpose,  we  consider  a Liapunov  function  given  by 

A rn 

V(yj^)  = Yj^Dyj^,  where  D is  a positive  definite  matrix.  Then  the 
forward  difference  AV(yj^)  = ~ given  by 

/\  T -'1’ 

AV(yj^)  = -yj^Eyj^,  where  -E  = A DA  - D.  If  E is  a positive 
definite  matrix,  then  the  solutions  will  be  asymptotically  stable. 
On  the  other  hand,  if  we  have  asymptotic  stability,  i.e.,  all 
eigenvalues  of  A have  modulus  strictly  less  than  one,  then  for 
every  positive  definite  matrix  E there  exists  a unique  positive 
definite  matrix  D that  satisfies  the  equation  a'^DA  - D = -E, 

A A 

Moreover,  if  for  some  real  number  p,  0 < p < 1,  all  the 

1 


eigenvalues  of  the  matrix 


A have  modulus  less  than  one. 
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then  to  every  positive  definite  matrix  E there  corresponds  a unique 
positive  definite  matrix  D that  satisfies  the  equation 

AdA  - (l-u)D  = -E.  (3.5) 

If  AV(yj^)  = -y^Eyj^  - - ~ 

The  matrix  A of  particular  interest  here  is  seen  to  be 
equivalent  to  a matrix  in  companion  form.  In  this  case,  for  the 
existence,  uniqueness  and  positive  definiteness  of  the  matrix  D of 
equation  (3.5)  it  suffices,  [12],  for  the  matrix  E to  be  positive 
semidefinite  and  not  identically  equal  to  zero.  This  remark 
allows  us  to  choose  particularly  simple  matrices  E;  our  purpose 
is  to  obtain  as  simple  a form  as  possible  for  this  matrix  D and 

A. 

for  this  it  is  convenient,  given  the  special  form  of  the  matrix  A, 
to  restrict  ourselves  to  certain  choices  of  E.  We  shall  represent 
the  unique  solution  D of  equation  (3.5)  in  the  form 

(3.6) 

where  a is  a scalar,  r = (r^,...,rj^)  and  r = (r^,...,rj^) 

tv  ^ 

are  N-dimensional  vectors,  Q = ('Ij^j)  ^^1^3  Q = (^j^j) 

N X N symmetric  matrices  and  Q = is  an  N x N matrix. 

The  substitution  of  (3.4)  and  (3.6)  into  eauation  (3.5) 
suggests  that  a particularly  simple  form  for  the  matrix  D can 
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be  obtained  if  the  matrix  E = E = (e^^)  is  chosen  to  have  zero 

entries  everywhere  except  the  elements  “ ^N+1  1' 

®1,2N+1  ®2N+1,1'  ®N+1,N+1'  ^N+1,2N+1'  ®2N+1,N+1  ®2N+1,2N+1' 

Moreover,  the  simplicity  of  the  structure  of  the  matrix  A allows 

for  the  vector  r"^  = r"^  = ( r , . . . , r,j)  and  for  the  matrices  0,0 

and  O to  be  related  q.  . = q.  . = a,  . for 

i 7^  j / i/j  = 1,...,N  and  for  = 0 for  i = 1,...,N  if  the 

following  equations  are  satisfied 


i,  j — 2,... ,N, 
%,1  ' i = 2 N, 


(l  + ib)ri+  (l-a)r,.^=0, 


.,N, 

i = 2,...,N, 


and  the  nonzero  elements  of  E are  selected  as 


®1,N+1  ®N+1,1  “ ^1,2N+1  ®2N+1,1 


= -[(^  cu  - ar^)  (1  + ^ b)  + (1  - 1 cr,  + 


2k  ~ 

, a b ~ 
c ^11 


c N 1 
(i-y)rj^]. 


(3.7a) 

(3.7b) 

(3.7c) 

(3.7d) 


(3.7e) 
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'N+1,N+1 


[(i  c)2a  -2ar^  i c + - (l-a)q^^j^], 


{3.7f) 


'2N+1,2N+1 


-[(^  c)^a  - 2ar^  ^ c + a^q^^  - (l-M)q^^j^],  (3.7g) 


1 


®2N+1,N+1  ^N+1,2N+1  ~ 2ar^  ^ c + a • (3.7h) 

From  equations  (3.7a-b)  it  follows  that 


,,  ^,-N+i 

'li.j  = 'l-**’  •’n.n- 

(i-j)' 

{3.8a) 

/ 1 ''  V i- 1 

<3i,i  = =11,1' 

i = 1 , . . . ,N 

(3.8b) 

‘!l,i  ' '^1,1' 

i = 1 , . . . , N 

(3.8c) 
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Now,  letting  V(Xj^)  = where  = (x^,  x^,  . . . , xj^)  , 

^ ® ^ i?  using  (3. 7a-d)  , 

(3.8a-c),  (3.9)  and  (3.10),  we  obtain  the  desired  Liapunov  function 
in  the  form 


V(Xj^)  = “(x°)^  + 

^ ^ N %,N-(i-l)  ^ ^^“^^‘^N,N-1^  ■ 

i-1  i 

, i ^k  , T , 

■ - a T I iJ  + 

N 

+ i Id  - 2 i P)-“+ha  t i + d a ^ b)q^^,) 

0 1 
1 X,  - X, 

,1  k k,  , 

• (cx  - a + 


N i-1 

+ 2 I I (cx^  - a 
i=2  j=l 


i-1  _ i 

^^k M ->  ,n-N+i^ 

T ^ ^ N ^N,N- (i-j) 


x^"^  - x^ 

- a ^ + 
^ k I ^ N N 


r ,,  ^ 'C.-N+ig/  i\2  ^ , 

1=1 

«,  i-1  i 

N I ^T_L  ■ - X,  - X,  -r 

, (1  - 2 i)-N-^.(JL_^)2  ^ , 

^=1  N 


(3.11) 


and  its  forward  difference,  divided  by 
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AV(x.  ) ^ „ 

— = [2a(b+y)  + ab  ^ + 3 + b^Y  + 


+ 2(1  + ib)(o-ab)((l  - 2 i + i 


N 2 H ^ N N ‘^'‘^N,N-1  ’’’ 


+ ^ '=’‘Jh,i'  + - <1  - 2 ^ >'>”81  (x“)2  + 


+ 5 - 2 i P)-""'(a  + i 1 


'1  * S ‘>>‘’n,i’  ->-  '»^  - <1  - 2 i p)"y)  + 


1 2[ac(l  + i b)  1 1(1  - 2 i p)-«"l(a  + i 1 


(1  1 i b)q^_^)((l  - f)  i c2 


+ bcY  - (1  - 2 - u)  + (1 


ac  (1  + - b)  ] + 


2 N 


+ 2(-Sa(l  lib)  - in  - 2 ^ ^ ^ '=>'In,n-1  b 


b (1  b ^ b)q„_ji;(l  - f,  i ac  - a2(l  + I b) 1 
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- abY  + (1  - 2 ^ M)  + (1  - 2 1 P ) a]  x 

+ 2[-aac  1 + 2a^c  1 [ (1  - 2 ^ ^ + 

N-1  _ 

+ (1  + ^ Wq„,ll  - ca9lx»  -ii ^ - 2V,V(X^), 

N 


N-1  N 
0 ""k  ■ ""k 


T 

Fy 


(3.12) 


where  the  off  diagonal  terms  of  the  Q matrix  must  be  related  by 


%,i-2  ~ %,i-l  _ a(i  _ 2 - ^N,N-(i-l)  %,N-(i-2)  ^ 


T 

N 


N 


T 

N 


^N,i-1  ^ 2"%,i-2 


+ (1  - 2 ;:ry) 


T , -N+i-1 


N '^N,N-(i-l)' 

i = 3, . . . ,N, 


(3.13) 


4.  A Liapunov  Functional  for  the  Scalar  Neutral 
Difference-Differential  Equation. ' 

The  results  of  Section  3 permit  us  to  obtain  a Liapunov 
functional  in  an  explicit  form  for  the  scalar  neutral  difference- 
differential  equation  (2.1).  For  this  purpose  consider  the  limiting 
process  described  by 


■JN.i  -5^ i i 
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AV(x,  ) 


. . -,  O ^ i-N+1  2iJ0 

lim  (1  - 2 - u)  ->  e 

N^-oo 


-T  < 0 < 0. 


These  formal  limits  applied  to  equations  (3.11),  (3.12)  and  (3.13) 

yield  the  equations 


V(x^)  = ax^(O)  + 


+ 2x  (0)  f (ae^^  ^®'^^^q(T+0)  + q(-0))(cx  (0)  - axMe))de  + 
^ J_T  ^ ^ 


0 fO 


-T  j 0 


(cx^(0)  - ax^(0)  )e^^  ^^’'’^^q('r+0-B) 


(4.1) 


(cx^(3)  - ax^(B))d3d0  + 


+ f°  e^^^6xj(0)d0  + e^^^Y  x'^(0)d0, 

J-T  ^ J-T  ^ 


V(x^)  = [2u(b+y)  + 3 + b Y + 


+ 2 (c-ab)  [aq (T) e + q(0)]]x^(0)  + 

+ [Yc^  - 3e"^^'^]xJ(-x)  + [Ya^  - Ye"^^  ^ ] x^  (-T ) + 

+ 2c[a-q(x)  - 2aq(0)  - a^q(T)e^^^  + bYl x^ (0) x^ (-x)  - 


2a[a-q(T)  - 2aq(0)  - a^q(T)e^^^  + by]  (0)  {-x ) - 


- 2acYx^(-x)x^(-x)  - 2yV(x^), 


(4.2) 


q'(-0)  - ae^^  (0  + x)  = - (b+2ii ) q (-0  ) - 


_^2vi  (*^  + ^)cq(0+T) 


(4.3) 


At  this  juncture,  it  is  convenient  to  introduce  the  notation 
defined  by  the  equations 


q (6  ) = p (x-0  ) e 


y (T-.0) 


and  to  let 


a = a + p(0)  + 2ae^^p(x)  + a^p(0)e^^^, 

B » ie"'  , 

9 = Ye>'’  , 


which  introduced  into  (4.1),  (4.2)  and  (4.3)  yields  the  equations 


V(x^)  = 5x^(0)  + e^^  e^^®3x^(0)d0 


+ e^^^Yx'^ (0)d0  + 

-T  ^ 
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+ (p(0)  + 2ae*^^p(T)  + a^p(0)e^^  Sx^(O)  + 

+ 2x. (0)  ^ (p (T+0 ) + ap(-9)e^^)* 

^ -T 

.(cx^(6)  - ax^{9))d0  + 

+ 2 f°(cx  (0)  - ax:(0))e'‘^^^‘^^^^’^p(3-0) 

■ -T j 6 ^ ^ 

• (cx^(O)  - ax^(0))d3dO,  (4.4) 

V(x^)  = [2(b+M)a  + e‘^^{B+b^Y)  + 

+ 2(b+u)p(0)  + 2e^^ap  ( t)  (b+p)  + 

+ 2e^^^a  (c+ay)  p (0)  + 2e*^  ^ (c+ay)  p ( r)  ] x^  (0)  + 

+ [Yc^e'^^  - 3e"‘^"]x^(-T)  + [e'^^a^Y  - e"‘^’^Y]  (-x)  + 

+ 2c  [cii+bYe'^^  ] x^  (0)  x^  (-T ) - 

- 2a  [oi+bYe^^  ] x^  (0)  x^  (- t ) 

- 2ace*^^Yx^  (-X ) x^  {-X ) - 

- 2mV{x^)  (4.5) 

and 
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p'(0)  - ae^^p'Ci-0)  = (b+M)p(0)  + 
+ (c+aM ) e'^  ^p  (T-0  ) . 


(4.6) 


Equations  (4.4),  (4.5)  and  (4.6)  have  been  obtained  by  formally 

taking  limits  on  the  Liapunov  function  for  the  discrete 
approximation  to  our  original  functional  differential  equation.  It 
is  easy  to  see  that  (4.4)  is  a well-defined  functional  on  and  a 

straightforward  although  laborious  computation  shows  that,  if 
(4.6)  is  satisfied,  (4.5)  represents  the  rate  of  change  of  the 
functional  (4.4)  along  the  solutions  of  the  scalar  neutral 
functional  differential  equation  with  initial  conditions  in  the 
domain  of  the  generator;  we  postpone  such  an  analysis,  and  use 
the  results  thus  far  obtained  as  motivation  for  the  method  of 
analysis  presented  in  the  next  section. 


5 . A Liapunov  Functional  for  the  Matrix  Neutral 
Difference-Differential  Equation. 

The  results  of  the  previous  section  suggest  a form  for  a 
Liapunov  functional  for  the  matrix  neutral  difference-differential 
equation  (2.1) -(2. 2).  For  this  purpose,  on  the  space  consider 

the  real  symmetric  quadratic  form 


V(iJ>)  = ‘p'^(0)M(p  (0)  + e^^ 


rO 


-t 


)'^(e)Re^^‘®<}>  (Q)d0  + 


+ e 


4)'^'  (0)Se^^^<|)  ' (0)d0  + 
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+ ())'^(0)  [P(0)  + e^^A'^p'^(T)  + e^‘’P(i)A+  ‘ a'^P  { 0)  A]  0 ( 0 ) + 

+ 2<j)'^(0)  I (P  (ui+T)  + e^^A'^p'^(-u)  ) • (C(|>  (a)  - A())'(a))du  + 


)'^(a)c'^  - <1)'^'  (a)  a'^)  P (6-01)  e*^*  • 


• (C4)  (B)  - A<j) ' (B)  )dBdu, 


(5.1) 


where  p is  a real  number,  M,R,S  are  constant  n x n real 
positive  definite  matrices,  and  P(u)  , 0 < cx  < t,  is  a continuously 
differentiable  matrix;  we  assume  that  F(^)  is  a solution  of  the 
initial  value  problem  for  the  functional  differential  equation 

pt(a)  - e'^^A’^p'^' (T-a)  = (B+pI)P(a)  + 

11 T T mm 

+ e^  (C  + yA  )P^(T-a),  o < ct  < i , 

P(0)  = P'^(O)  = Pq, 


(5.2) 

(5.3) 


where  Pq  is  an  arbitrary  symmetric  matrix. 

Consider  the  Frechet  differentiable  functional  (5.1)  evaluated 
along  a solution  of  (2.1)-(2.2)  with  initial  conditions  in  ; 

this  yields  a function  of  time,  denoted  by  V(t)  = V(x^) , which  is 


differentiable  along  such  solution ; after  a straightforward  but 
lengthy  computation,  and  through  use  of  (5. 2) -(5. 3),  we  obtain  that 
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V(t)  = V(x^)  = V(x^)  = -2mV{x^)  + 


T T VI C T 

+ x^(0)  I (B  +MI)M  + M{B+MI)  + 2e^^(R+B  SB)  + 


+ (b'^+mI)P(0)  + P(0)  (B+yl)  + (c'^+Ma'^)p'^(T)  + 


+ e^^P  ( t)  (C+UA)  + '■A'^p'^(t)  (B+mD  + 


+ (b'^+mI)P(t)  A + e^^^A'^P(O)  (C+MA)  + 


2u  I T T 

+ e ^ (C  +yA  )P(0)A]x^(0)  - 


- h(x^(0)  ,x^(-t)  ,x^(-T)  ) = U(x^), 


(5.4) 


where 


h(x^(0),  x^(-T),  \(-T))  = [x^(0),  -e  ^'^x^(-T),  -e"^^x^(-T)] 


R + b'^SB 

c'^(M+e^^SB) 

R - 

T U T 

A (M+e^  SB) 

e^'-^A'^SC 

2yT  T 


X^(0) 


-e  ^"^x^C-x) 


S - e^^^A'^SA  -e"^^x^(-T) 


(5.5) . 

This  computation  is  valid  only  on  solutions  with  initial  data' 


on  3I{q/);  however  a direct  application  of  Theorem  3.9  of  [16]  yields 
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that  for  any  solution  with  initial  condition  on  we  have 

that  for  V(x.)  = lim  r < U(x  ). 

^ A\0  ^ 

The  objective  of  this  paper  is  to  show  that  it  is  possible 

to  give  an  estimate  of  the  rate  of  growth  or  decay  of  the  solutions 

of  the  matrix  neutral  difference-differential  equation  (2.1) -(2. 2), 

using  a functional  of  the  form  (5.1)  and  its  derivative  (5.4). 

For  this  purpose,  let  Y = max  { Re  | A G o (_a/)  } . Given  e > 0 

and  -p  = Y + 2e  we  wish  to  show  there  exist  matrices  M,R,S  and 

a differentiable  matrix  P(a)  satisfying  (5. 2) -(5. 3)  so  that  for 

the  functionals  V(<}>)  and  U(<i))  given  by  (5.1)  and  (5.4)  with 

these  matrices  we  have  that 

cj  |4>M^  < V(<1))  < c^\\<P\\^  (5.6) 

and 

V(<J))  < -2MV(<i))  , (5.7) 

for  some  positive  constants  c^,C2-  If  this  is  possible,  a norm  is 

induced  by  the  square  root  of  the  Liapunov  functional  (5.1), 
which  we  denote  by  = {V (<(>)} ^'^^ . Relationships  (5.6) 

and  (5.7)  show  that  the  norms  ^ ^ ‘ ^ equivalent 

on  iif,  and  that 
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or 


I |x^| 


nl/2 


-u  t 


(5.9) 


The  estimates  (5.8)  and  (5.9)  are  precisely  those  stated  in  (2.4). 


Note  that  the  norm  | | * | I best  possible  one  in  the  sense 


that  it  yields  (2.4)  with  K = 1.  Also,  if  'i  < 0,  it  follows 
from  (5.8)  or  (5.9)  that  the  Liapunov  functional  (5.1)  proves  uniform 
exponential  asymptotic  stability  for  the  solutions  of  (2.1). 


In  this  manner,  we  have  the  following. 


Theorem  1:  Consider  the  matrix  neutral  difference-differential 


equation  with  one  delay 


x(t)  + Ax(t-T)  = Bx(t)  + Cx(t-T) 


1 


and  the  Liapunov  functional  given  by  equation  (5.1).  Let 

Y = max{Re  X|det[X(I+Ae  ^^)  - B - Ce  = 0} 

and  ^ > 0.  If  there  exist  constant  positive  definite  matrices  M,R,S 
and  a differentiable  matrix  P (oc ) , 0 < a < t , that  satisfies  (5.2) 
with  P(0)  = P(O)'^,  such  that  the  inequalities  (5.6)  and  (5.7)  hold 
with  Positive  and  with  -u  = Y + 2^;  then  the  solutions  of  the 

difference-differential  equation  satisfy  the  exponential  bound  (5.9) 
and  the  equation  is  exponentially  asymptotically  stable  if  u > 0. 

In  the  next  two  sections  we  show  that  an  appropriate 
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differentiable  matrix  P(ot)  as  required  by  the  above  theorem  can 
always  be  chosen,  and  we  analyze  its  structure.  In  a further 
section  we  determine  appropriate  matrices  M,R  and  S so  that  the 
conditions  for  the  above  theorem  are  always  satisfied.  We  thus 
demonstrate  the  existence  of  a Liapunov  functional  of  the  stated  form 
for  the  neutral  functional  differential  equation. 

At  this  juncture,  we  remark,  for  repeated  later  use,  that  for 
the  particular  above  choice  of  -y  = y + 2e,  the  spectral  radius  of 
e^^A  is  always  less  than  one.  Indeed,  letting  x(t)  = e ^^y(t), 
t > 0,  our  original  functional  differential  equation  (2.1)  becomes 

y(t)  + e^^Ay(t-i}  = (B+Ml)y(t)  + (C+uA)  e^  V ) / t > 0, 

for  which  the  solution  operator  is  given  by  T(t)e^^,  where  T(t) 
is  the  solution  operator  of  equation  (2.1).  Now,  using  equation 
(2.4)  and  -y  = Y + 2^,  we  note  that  | |T(t)e^^|  | < Ke 

e > 0,  which  implies  the  uniform  asymptotic  stability  of  the  solutions 
of  the  above  equation.  But,  from  [5] , this  implies  that  the 
solution  to  the  difference  equation  y(t)  + e^  Ay(t-i)  =0  is  also 
uniformly  asymptotically  stable.  But  this  can  be  the  case  if  and 
only  if  the  spectral  radius  of  e^^A  is  strictly  less  than  one. 

6.  The  Functional  Differential  Equation  for  P(a). 

In  this  section  we  consider  the  functional  differential 


equation 


P'  (a)  - e^^A'^p"^(x-ai)  = (B'^+pI)P(a)  + 
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+ (c'^+Ua'^)  p'^(x-a)  , -CO  < a < oo, 


(6.1) 


with  the  initial  condition 


P(2)  = K, 


(6.2) 


where  K is  an  arbitrary  n n matrix.  V7itliout  any  loss  of 
generality,  let  p = 0 in  equation  (6.1),  thus  oljtaining 


P'(ot)  - a'^p"^(t-u)  = B'^P(oi)  + 


T T 

+ C P (x-a)  , _oo  < a < 


(6.1' ) 


For  the  particular  case  A = 0,  namely 


P'(a)  = B'^P(ot)  + c'^p'^(x-oi)  , -oo  < a < «>  , 


(6.1") 


it  has  been  shown  in  [1]  that  the  solutions  of  such  an  equation 

with  initial  conditions  given  by  (6.2)  exist  and  are  unique  and 

that  the  linear  vector  space  of  all  solutions  of  such  an  equation 
2 

has  dimension  n ; moreover  an  algebraic  representation  of  these 
solutions  was  presented.  The  results  for  the  more  general  equation 
(6.1')  are  analogous  to  those  presented  in  [1]  for  equation  (6.1") 
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I 


i 

i 


and  we  only  sketch  them  here.  We  assume  in  the  sequel  that  all 
the  eigenvalues  of  A lie  strictly  inside  the  unit  circle. 

Defining  the  matrix  F(a)  = p'^(c-ot),  the  functional  differ- 
ential equation  (6.1')  with  the  initial  condition  (6.2)  reduces  to 
the  system  of  ordinary  differential  equations 


P'  (a)  + a'^F'  (a)  = b'^P(u)  + c'^F(a) 

P'(a)A  + F'  (a)  = -P(u)C  - F(u)B 


(6.3) 


or 


P'  (a)  - A'’^P'(a)A=  B'^P(a)  + A'^P(a)C  + c'^F(u)  + A'^F(a)B 
F'  (a)  - a'^F'  (a)A  = -P(a)c  - B'^P(a)A  - F(a)B  - c'^F(a)A 


(6.3'  ) 


with  initial  conditions 


P(j)  = K,F(j) 


Consider  the  notation 


= K 


(6.4) 


P(a)  = (p^j(a))  = 


PT*(a) 


where  P^*(a)  and  p* ^ (a)  are,  respectively,  the  i^^  row  and  the 
column  of  P(a).  Then,  equations  (6.3')  and  (6.4)  can  be  rewritten. 


through  the  use  of  the  Kronecker  (or  direct)  product  of  two  matrices 
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2 

[10]  as  the  2n  system  of  ordinary  differential  equations 


I - A g»  A 

0 

d 

P (w) 

0 

I - A (gi  A 

dot' 

f (cx) 

B(8iI+A®C  C!X»I  + A<X)B 

-I(8)C-B«>A  -I»B-C(»A 

L_ 

with  initial  conditions 

p(i)  = f(“>  = ikli «;•';) 

Using  the  fact  that  1 is  not  an  eigenalue  of  A,  then 
I - A is  invertible  and  so  is  I - A ® A.  Thus  the  system  (6.5) 
can  be  rewritten  as 

dot 

(I  - A (g)  A)”^- 

• (B  0 I + A (8>  C) 

(I  - A ®A)“^- 

• (-1  (g)  C - B ig>  A) 


(I  - A <g)  A) 

• (C  ig)  I + A <g>  B) 

(I  - A A)"^' 

• (-1  (gi  B - C <gi  A) 


p{ot) 


f(a) 


(6.5’) 


p(“) 
f (ct) 


p (ot) 

f (u) 


(6.5) 


Consider  now  the  uniqueness  of  solutions  of  ( 6 . 1 ' ) - ( 6 . 2 ) . If 

k 


I 


t 


i 

1 


1 

I 


ti 


!' 

I 

I 

I 


this  initial  value  problem  has  a differentiable  solution  P(ct) 
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defined  on  -t”  < a < <»,  it  follows  that  the  pair  of  matrices 
P(a),F(ct)  with  F(a)  = p’^(T-a),  will  satisfy  (6.3),  (6.4) 

and  then  the  pair  of  vectors  p(«),f(a)  as  defined  above,  will 
satisfy  (6.5') > (6.6).  Because  of  the  linearity  of  all  the 

equations,  and  the  uniqueness  of  the  solutions  of  (6.5'),  (6.6), 
it  follows  that  the  solution  P(oi)  of  (6.1'),  (6.2)  is  unique. 

Consider  now  the  question  of  existence  of  solutions  of 
(6 . 1 ' ) - (6 . 2 ) . The  initial  value  problem  (6.5')-(6.6)  has  a unique 
solution  (p(a),f(a))  defined  on  -oo  < a < <».  This  implies  the 
existence  of  a unique  pair  of  differentiable  matrices  P(<^),F(u) 
defined  on  -oo  < a < » and  satisfying  (6.3),  (6.').  These  equations 

can  be  rewritten  as 

^ P(a)  + ^ F(a)  = b'^P(u)  + c"F(a) 

^ f'^(t-u)  + a”^  ^ p'^(T-a)  = b'^F(t-u)  + c'^P(T-a), 

with  the  initial  condition 

t T l 

P(^)  = K = F^  (-)  . 

Then  it  follows,  from  the  uniqueness  of  the  solutions,  that 
T 

F(a)  = P (T-a)  for  < u < «. 

Moreover,  it  is  well  known,  [2,4],  that  the  system  (6.5')  has 

2 

2n  linearly  independent  solutions  which  can  be  represented  in  the 
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following  fashion.  Let  £ < 2n  , be  the  distinct 

eigenvalues  of  the  matrix  in  (6.5' )»  that  is,  solutions  of  the 
xdeterminental  equation 


det 


(I  - A g)  A)  X - (D  g I + A g C) 


I g C + B g A 


-C  g I - A g E 


(I  - A g A)X  + (I  g B (}>  C g A) 


= 0;  (6.7) 


each  X.,  j = 1,...,?,,  v/ith  algebraic  multiplicity  m.  and 

^ r ^ r 2 ^ 

geometric  multiplicities  n.,  5!n.=m.,  Im.  =2n.  Then,  the 

3 r=l  3 3 j J 

2 

2n  linearly  independent  solutions  of  (6.5')  iire  given  by 


<p?  (a)  = e 


T q-i 

_2j 

■ "j 


X (a-i)  a - i)'!-  ^ 

L nr _ , 


(6.8) 


2 

where  q = l,...,n^,  r = l,...,s,  I n.  = m. , I m.  = 2n  , and  the 


3 


r=l  ^ ^ 3 ^ 


2n'^  linearly  independent  eigenvectors  and  generalized  eigenvectors 
are  given  by 

[X  l-H]e^  = -e^"^  , e°  = 0, 

where  H is  the  2n^  x 2n^  matrix  in  (6.5').  Changing 
notation  and  returning  from  the  vector  to  the  matrix  form,  we  see 
that  the  2n^  linearly  independent  solutions  of  (6.3)  are  given  by 


3fr 


(a) 


V?  («) 

. 3,r'  'j 


= e 


X . (C  - f 


q 

I 

i=l 


(q-i) ! 


3»r 


L 3»r. 


(6.9! 
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for  q = r = j = 1 = m., 

3 n=l  ^ J 

2 

F m . = 2n  , where  the  generalized  eigenmatrix  pair  (l1  ,M^  ) 

1 3,r  j,r 

associated  with  the  eigenvalue  satisfies  the  equations 

i i T T i i T i-1  T i-i 

(X  ^ - AlJ  ^(cSxa')  - CmJ  - Af1^  b'  = AL.  " L . ^ , 

(6.10) 

+ BL^  a"^  + (A.I+b"^)  + (C-Xa)M^  a'^  = ~ 

for  i = l,...,nT,  r = 1,  . . . ,s,  L*?  =0,  M*?  ^ = 0.  The  structure 

J 3 J » = 

of  these  equations  is  very  special;  in  fact,  if  they  are  multiplied 
by  -1,  transposed,  and  written  in  reverse  order,  we  have 


.T  .T 


(-X  I-B)M^"  - AM^^  (c'^-aa^)  - CL^  - AL^  B^  = -AM^  ^ a"^  + J;  , 

D D,r  j,r  :,r  j,r  j,r  j,r  ^ 

■jf  T T i”^  T i"^  T i-l'^  T i-l'^ 

m!-  C + BMi-  A^  + L (-X.I+B^)  + (C+XA)L  a"  = -AL.  ^ A + L . ^ , 
D,r  j,r  :,r  j j,r  j,r  j,r 


(6.11) 


qT  qT 

for  i = l,...,n^,  r = 1, . . . ,s,  L ^ = M ^ = 0,  This  result  shows 
3 3 3 

that  if  Xj  is  a solution  of  (6.7),  -X^  will  also  be  a solution; 

moreover,  X^  and  -X^  have  the  same  geometric  multiplicity 

and  the  same  algebraic  multiplicity.  Thus,  the  distinct  eigenvalues 

always  appear  in  pairs  (Xj,-Xj).  An  examination  of  equations 

(6.10)  and  (6.11)  shows  that  if  the  generalized  eigenmatrix  pair 

corresponding  to  X , are  (L^  ' then  the  generalized  eigen- 

3 J»rj,r  ,p  '1''^ 

matrix  pair  corresponding  to  -X^  will  be  ((-1)^  ^j,r'  ^j,r^‘ 

These  remar)<s  show  that  if  the  solution  (6.8)  corresponding  to  X^  is 
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added  to  the  solution  (6.8)  corresponding  to  multiplied  by 

(-1)'^'''^,  we  have  n^  linearly  independent  solutions  of  (6.3)  given  by 


5?  (u) 

3»s 


TT?  (a ) 


X . (u 


i=l  i 

M . 


-X  . (a  - 


7)  q (^ 

2 I — 

i=l 


_ l\q-i 


.T 


(q-i) ! 


.T 
3 »r 


£ T 

which  satisfy  the  condition  E?  „ (y)  = This  condition 

j,rz 


is  precisely  condition  (6.2);  it  therefore  follov/s  that 


q (u  - il''  '■  i 

=lr'“>  ■ (q-i)!  “ q,r  " 


-X . (a  - -)  .T 
4-  (-D^+^e  3 2 1 

3 f r 


(6.12) 


for  q = l,...,n^,  I n^  = m . , m.  = 2n^,  are  n^  linearly  in- 


^ J r=l  ^ 2j  3 

dependent  solutions  of  (5.10').  Hence,  we  have  shown  that 


equation  (6.1')  has  n linearly  independent  solutions  given  by 


equation  (6.12),  where  the  generalized  eigenmatrix  pairs  (L^  ,M^  ) 

3 , s 3 f s 


satisfy  equation  (6.10)  for  one  of  tlie  elements  of  the  pair 


(Xj,-Xj),  each  of  which  is  a solution  of  equation  (6.7). 
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Making  use  of  the  above  results,  it  is  clear  that  the  functional 
differential  equation  (5.2)  with  initial  condition  (5.3)  has  a unique 
solution  P(o()  defined  for  -«  < a < <». 

7.  The  Matrix  Function  P(a)  and  Its  Properties. 

In  this  section  we  show  that  the  solution  P(a)  of  (6.1) 
has  an  integral  representation  which  is  particularly  useful  for  our 
purposes.  To  this  end  consider  the  matrix  function  P(^)  defined  by 

-OO 

P(a)  = Y'^(3)e^‘^WY(3-^^)e'^  ct  > 0 , (7.1) 

0 

where  W is  an  arbitrary  symmetric  matrix,  Y(3)  is  the  solution 
of  equation  (2.6)  and  y = -y  - 2£,  fc.  > 0.  This  integral 

always  converges  since,  for  every  e > 0 » we  have  | |y(3) | | < 

Ke(Y+£)3  jqj,  K > 1. 

From  (7.1)  it  immediately  follows  that 

P(0)  = P^(0) 

and  that 

P(a)  = Y'^(3)e^^«(3-‘t)e*' ^*^'^^d3  = Y'^(uH)e^  ^^‘''^\'JY(^)e''^d5 

Ja  0 

= p'^(-a),  _co  < a < 00.  (7.3) 

We  now  show  that  P(a)  satisfies  the  equation 

p'  (a)  - '^(  t-ot)  = (B'^+Ml)P(a)  + 

+ (C^+UA'^)e^ 'p'^(x-«)  , 0 < « < T. 


(7.4) 
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Indeed,  using  equation  (2.6)  in  the  definition  of  P(u),  we  have 


P(a)  = 


(I+A)dg  - 


0 

■“ 

0 

r 

0 


y'^  ( 6 ) ^ ^ Y ( - 1; ) Ad  3 


■& 
a 
'3-t 
a 


Y'^(3)e^^fe^^^~^^Y(3-C)Bd5d3  + 

Y^  (3)  e^  ) Y ( 3-C-^ ) CdCd3 . 


The  order  of  integration  in  the  last  two  terms  can  be  interchanged, 
yielding 


P(a)  = e 


-ya 


- e 


Y'^(3)e^^We^*^(I+A)d3  - 


Y'’’(3)e^^Wd^^^'^”^^Y(3-«-'t)Ad3  + 


r 

a 


Y'^(3)e^^^^^^"^^Y(3-C)d3 


y"^  ( 3)  e^  ^ Y ( 3-S-t ) d3 


or 


p(a)  = [ Y’^(3)e^‘^We'‘^(I+A)d3  - e^^P(a+T)A  + 


+ e“^“  I [P(5)e*'^B  + P(5  + x)e^‘ 

a 


(7.5) 
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Every  term  in  (7.5)  is  an  absolutely  continuous  function;  hence 
differentiation  of  (7.5)  gives 

i-  P(a)  = P(u+T)]A  - P(a)B  - e'^^P(a+T)C  - 

doc  doc 

-00 

- Y'^(e)e^^We^‘^(I+A)d3  + 

Jo 

+ e"^"^  [P(C)e^^B  + P(^+i)e^^^‘^^^C]dC) 

• a 

or 

i_  p(a)  = P{a+'C)]A  - P(a)B  - e^'^P(a+T)C  - 

doc  doc 

- M(p(a)  + e^ '■p  (ot+ t)  A) 


or 

— P(u)  + P(a+T)]A  = -P{oi)  (B+yl)  - 

doc  dec 

- P (a+T)  (c+yA) e^^ 


p'  (a)  + e^^^P’  (w+'C)A  = -P{oi.)  (B+yl)  - P (u+i: ) (C+yA)  e^  ^ . 
We  have,  therefore,  that  P(ot)  satisfies  the  equation 


P'  (-a)  + e^^P'  (-a+'c)A  = -P(->x)  (B+yl)  - P(-<-t+T)  (C+yA)e^^; 
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but  using  property  (7.3),  we  obtain  that 


P'  (-a)  = - 


^ P(-a)  = 
da 


d ~T 


~ ' T 

-P  (a) 


Hence 


-p"^(a)  + e^^P'  (-a+t)A  = -p'^(a)  (B+til) 


- p(-u+T) (c+MA)e 


U 't 


or 


P'(«)  - e^^A'^p"^(x-a)A  = (B'^+Ml)P(a>  + 

UT  T T ~T 
+ e^  (C  +MA  )P  (T-u) . 

Because  of  the  uniqueness  of  the  solutions  of  (5. 2) -(5. 3), 
it  is  seen  that  P(a)  defined  by  (7.1)  is  the  unique  solution 
of  (5. 2) -(5. 3)  with  the  initial  condition  prescribed  by  (7.2). 

We  now  proceed  to  show  that  the  following  relationship 
is  satisfied: 


Z = (B  +MI)P(0)  + P(0)(B+uI)  + 

+ e^^  (c'^+MA^)  P^(x)  + e*' 'p  ( I ) (C+mA)  + 

+ e^^A’^p'^(T)  (B+nD  + e^^  (B'^+ViI)P(T)  A + 

+ e^^^A'^P(O)  (C+yA)  + e^^ ^ (c'^+ua'^)  P ( 0)  A = -W. 


(7.6) 
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Indeed,  using  P(0)  and  P'T)  given  by  (7.1)  and  the  relation 


^P(O)  = Y'^(3-T)e^^WY(3-'t)e^^dB, 

n 


it  is  easy  to  see  Z can  be  represented  as 


mOO 

Z = [b'^Y(6)  + cV(e--'C)  ]G^^We^*^'[Y(B)  + Y(3-OA]d3  + 

0 


+ f [y'^(B)  + A'^Y'^(3-'t)  ]e‘^^'}e^^[Y{3)B  + Y(3-T)C]d3  + 


.00 

+ 2y  [y'^(3)  + A'^Y'^(3-t)  ]e^*^We^^[Y(3)  + Y(3-’t)A]d3. 
n 


Making  use  of  (2.3),  we  have  that 


.00 

= ^ [[y'^(3)  + a'^y'^(3-t)  ]e^®VVe'^^[Y(3)  + Y ( 3-1: ) A]  ] d3 , 


from  where  it  follows  that  Z = -W. 

We  remark  that  to  each  constant  symmetric  positive  definite 

matrix  W there  corresponds  a matrix  function  P(ct)  given  by 

(7.1).  This  matrix  P(ct)  is  the  unique  solution  of  (5.2)-(5.3) 

~ T 

with  the  initial  condition  prescribed  by  P(0)  = P(0)  = 

.00 

Y'^(S)e^^WY(3)e^^d3  which  is  necessarily  positive  definite. 

0 »p 

Conversely,  to  each  P(0)  = P (0)  = Pq,  equations  (5. 2) -(5.3)  have 
a unique  differentiable  solution  P(oi)  which,  through  (7.6)  defines  a 
unique  symmetric  matrix  W with  the  property  that  for  this  W, 
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(7.1)  gives  the  unique  solution  for  (5.2)-(5.3)  with  the  pre- 
scribed initial  condition.  It  follows  that  the  map  W ->■  P(0), 
defined  by 


P(0)  = f Y'^{3)e‘' V/(3)e^'^d3, 
J 0 


as  a map  on  the  space  of  n x n symmetric  matrices  is  one-to-one  , 
onto,  and  it  maps  positive  definite  matrices  W into  positive 
definite  matrices  P(0). 

8 . Structure  of  the  Liapunov  Functional  and  Its  Derivative. 

The  characterization  given  in  the  previous  sections  for  the 
matrix  function  P(ct)  permit  us  to  put  into  evidence  the  particular 
structure  of  the  Liapunov  functional  (5.1)  and  its  derivative  (5.4). 
In  fact,  the  substitution  of  the  representation  (7.1)  for  P(x)  into 
the  Liapunov  functional  (5.1),  yields 


V((I))  = 4>'^(0)M(f  (0)  + e*^^ 


fO 


(j)'^(0)Re^^°(})  (0)de  + 


-T 


+ e^^  r (t>  "^(‘3)Se^^®<t)  ' (0)d0  + 

-T 


+ I |[Y(3)  + Y(3-'c)  A]4>  (0)  + 

)T  2uB 

Y ( 3-a- i ) [C})  (a)  - Al)'(a)]dix|  We  . 
-T  ^ 

{ [Y(3)  + Y(3-OA]<))  (0)  + 


Y(3-«-t)  [C4)  (tt)  -A(p'(a)]da| 
-I 


d3. 


But  this  functional 


using  (2,5),  can  be  represented  as 


V(x^)  = x'^(0)Mx^(0)  + e''*^ 


rO 


-T 


x'^  (0  ) Re  ^^x^  (0  ) dO 


+ e 


Irp  2u0  ’ 

X ^ (0)Se  ^ X {0)de  + 

. T ^ t 


(8.1) 


where  x,^  d(0)  is  the  solution  of  (2.1) -(2. 2)  for  B > 0 \'/ith 

t+P  ~ 

initial  condition  x^(0),  < 0 < 0. 

In  an  analogous  manner,  using  the  above  notation,  equation  (5.4) 
can  be  rewritten  as 


V(x^)  = -2mV(x^)  + 


+ x'^(O)  [-W  + (A'^+Ul)ri  + M(A+yI)  + 2e^'^  (R+b'^SB)  ]x^(0) 


- e^^[x'^(0),  - e^^x'^(-i),  -e  ^^x^(-t:)] 


R + B SB 


T Li  T 
C [M+e^  SB) 


[M+e‘'  ‘b’’^S]C 


2u  T T 

R - c C SC 


T , U t , 2u  I T 
-A  [M+e  SB)  e A SC 


U t T 

-[n+e  B S]A 


2m  T T 
e ^ C SA 


S - e^t'^'^SA 


x^(0) 

-e'^\^(-T) 


-e  x^(-T) 


(8.2) 
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9 . Bounds  for  the  Liapunov  Functional  and  Its  Derivative. 

The  structure  of  the  Liapunov  functional  (5.1)  and  its 
derivative  (5.4)  developed  in  the  previous  sections  allows  us  to 
establish  that  it  is,  always  possible  to  choose  positive  definite 
matrices  M,R,S  and  W (and  therefore  P(oi))  so  that  bounds  of 
the  form  (5.6)  and  (5.7)  always  hold.  For  this  purpose,  first  we 
will  show  that  there  exist  positive  constants  that 

Cil  |4>  I ^ 1 V(9)  < C2I  |9  I 1^,  (^'1) 

for  all  (f)  G ^(_a/)  . 

Denote  bv  A . (H)  and  A (M)  the  minimum  and  maximum 
min  max 

eigenvalues  of  the  positive  definite  matrix  M,  the  corresponding 
notation  to  be  used  for  the  positive  definite  matrices  R,S  and  W. 

From  (8.1)  it  iirunediately  follows  that 


min(Amin(M),e-|^I^A^.^(R),e-i^l^^,^(S)) 


Using  (2.4),  it  is  easily  seen  that 


11^  : v(*: 


Hgio)  1 1 ; 11^1 1 


(A  +£)  3i  I 
< Ke  ' X 


O' 


Thus,  from  (8.1)  it  follows  that 


1 
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V (9  ) < [max  (A  (M)  ,e  ' 
max 


•A  __  (R)  ^A__.(S)  ) + 


max 


max 


, ^max  „2 , I I ^ I I 2 

+ 2 e K ] I I <*>  I I ^ 


The  desired  values  of  will  therefore  be  given  by 


^1  = 


|u(  ^ 


and 


c„  = 


max  (A  . (M),e^^^''A  (R),e^'^l'^ 

min  ' max  ' 


A (VJ) 

X (S))  + -i'lll k2, 

max  2 e 


Secondly,  we  proceed  to  show  that  enuation  (5.7)  always  holds. 
Equation  (8.2)  shows  that  we  only  need  to  choose  the  positive 
definition  matrices  W,R,S  and  M in  such  a manner  that  the  function 


h^(x^(0)  ,x^(-T)  ,x^(-i) ) = [x|^(0)  ,-e  ^\^(-0  ,-e  ’'‘'■x'^(-'r)  ] • 


T U 1 T U T U T T U T 

\V-(A  -HJI)M-M(A-HiI)  [M+e  B S]Ce  -[Wfe  B S]Ae 

T U'C 
-(R+B  SB)e 


T,  ^ 

C [^^  SB]e 


T,  , yi: 

-A  [Wte  SBle*^ 


[R-e^^‘  V^SC]  e^^  ^ e^*^*  SA 


3ut  T 
e A SC 


, 2111  T ,111 

[S-e  A SAle 


x^(0) 


-111  T, 

-c  x^(-0 

-e  x,_(-i) 


is  nohnegative.  We  note  that  this  form  can  be  rewritten  as 
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hj^{x^(0)  ,x^(-x)  ,x^(-T) ) = 


= e^’^[x^(0)  (rtt€‘^Vs),-e  ^'•x^(-x)c'^,-e  ^^x^(-:)a'^] 


+ X^(0)  [W  - (a'^+UDM  - M(A+MI)  - (R+G'^SB)e  - 

- a + ~ s“V  ^^^)(rB«^^SB)]x  (0)  + 

1 ^ 

+ e“^^x^{-T)  [R  - c'^(Se^'^^  + I + |-  Se^^^)C]x^(-T) 

+ e ^^x^(-x)  [S  - (l+£^)e^^^A'^SA]x^(-x)  , (9.2) 

where  is  an  arbitrary  positive  number. 

First  of  all,  note  that  since  the  spectral  radius  of  e^^A 
is  strictly  less  than  one,  then  for  > 0 sufficiently  small 

there  always  exist  a positive  definite  matrix  S such  that 

S - (1-t  e. ) ^a'^SA  > 0. 


'1 

4 

'! 
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R is  then  a positive  definite  matrix. 

The  matrix 

W - (A'^+m)M  - M(A+mI)  - {R+B'^SB)e^^  - 

- e^"^  (Mfe^^B’^S)  (I  + |-  s“^e”^^'^)  (W^SB) 

1 

will  certainly  be  positive  definite  bv  takinq  ri  = r and  w = k I 

W ’ 

> 0 sufficiently  large. 

Finally,  from  the  particular  form  of  the  first  term  in  (9.2) 
it  is  easily  seen  that  this  term  is  positive  semidef inite . 

In  this  manner,  it  is  seen  that  the  form  h^  (x^  ( 0 ) , x^  ) , x^.  {-l ) ) 
can  be  always  made  nonnegative. 

All  the  above  results  can  thus  be  summarized  in  the  form  of  a 

Theorem  2 ; Consider  the  matrix  neutral  difference-differential 
equation  with  one  delay 

x(t)  + Ax(t-x)  = Bx(t)  + Cx(t-T:) 

and  the  Liapunov  functional  given  by  equation  (5.1).  Let 

” X t " X T 

Y = max {Re  X|det[X(l+Ae  ) - B - Ce  ] = 0} 

and  e > 0.  Then  there  exist  constant  positive  definite  matrices 

M,R  and  S,  and  a differentiable  matrix  P(a),  0 < a < t, 

T ... 

with  P(0)  = P(0)  such  that  the  functional  V is  positive 


definite,  bounded  above,  and 


V < 2 (Y+t-)  V. 

Note  that,  if  Y < 0,  then  the  above  result  yields  ex- 
ponential asymptotic  stability. 
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